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Abstract 

We study clean superconducting quantum dots (SQD) and also site dependent Josephson cou- 
plings, on site charging energies and the intersite interactions in presence of gate voltage. We 
predict the existence of different fractionally quantized Cooper pair stair case with many inter- 
esting physical properties. The appearance of stair case is not only due to the Coulomb blocked 
phenomena but also for the site dependent Josephson couplings. We also explain physically the 
absence of other fractionally quantized Cooper pair stair case. The physics of fractionally quan- 
tized Cooper pair stair case has close resemblance with the fractionally quantized magnetization 
plateau of a spin chain system under a magnetic field. 

Keyword's: Superconducting quantum dots, Fractionally quantized plateau, Abelian bosoniza- 
tion and Renormalization group study 
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1. INTRODUCTION 



In last few decades mesoscopic physics of nanoscale superconducting systems are 
revealing many interesting properties. There has been intense research aimed at developing 
the superconducting flux based digital electronics and computers . The single Cooper 
pair box (SCB) and single Cooper pair transistors have developed experimentally and used 
to demonstrate the quantization of Cooper pair on a small superconducting island jjj, 
which is the foundation of charge qubit [J]. 

Here we mention briefly the mechanism of Cooper pair quantization for SCB (= Joseph- 
son junction circuit consisting of a small superconducting island connected via a Josephson 
tunnel junction to a large superconducting island). The SCB was first experimentally real- 
ized by Lafarge et al. |5| observing the Coulomb stair case with the step of 2e. The evidence 
of Coulomb stair case has also predicted in Ref. 6|. Realization of the first charge qubit by 
manipulation of SCB and the observation of Rabi oscillation was done by Nakamura et al [7]. 
Our operating point is in the charge regime, i.e. E c » Ej. We consider a finite numbers 
of Cooper pairs (n) in the mesoscopic island. The eigen state equation of SCB is 

E c {h — n g ) 2 \n >= E n \n > (1) 

n = 0, 1, 2 corresponds to the charge state with energy spectrum E n = E c (n — n g ) 2 . n g is 
the gate voltage induced charge in SQD. From Eq. 1, it is clear to us for a specific value of 
gate voltage (i.e n g = 1/2) the charge state |0 > and |1 > become degenerate. Switching on 
a small Josephson coupling will then lift the degeneracy and forming a two level system. So 
the system reduced to the qubit system. 

H SC B = -l/2(B z a z + B x a x ), (2) 

where B z = E c (l — 2n g ) and B x = Ej. The qubit level energies are given by the equation 



Ei t 2 = ± (l/2)yi? c (1 — 2n g ) + Ej . At the charge degeneracy point, n g = 1/2, at this 
point only the off diagonal part contribute and the energy levels are separated by Ej and 
the qubit eigen states \E\ >= |0 > — 11 > and \E 2 >= |0 > +|1 >. For these states, the 
average charge on the island is zero, while it changes to ±2e from the degeneracy point, 
where the qubit eigen state approach the pure charge state. 
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It is clear to us from the analyses of previous studies that the Cooper pair stair case 
(charge quantization state of SCB) only occurs due to presence of Coulomb blocked phe- 
nomena. Here we consider a model systems consists of array of clean superconducting 
quantum dots (SQD) and also a SQD systems with site dependent Josephson couplings, 
on-site charging energies and the intersite interactions. We will study different fractionally 
quantized Cooper pair stair case. It is clear to us from Eq. 1 that the basic physics of SCB 
can be understood by the spin chain model Hamiltonian under a magnetic field. We will 
see from the analysis of our model that the factionally quantized Cooper pair stair case is 
nothing but the fractionally quantized magnetization plateau state of spin chain system. 

2. COOPER PAIR STAIR CASE AND PHYSICAL ANALYSIS OF COOPER OF 
CLEAN SUPERCONDUCTING QUANTUM DOTS 

The model Hamiltonian of SQD system has different parts, 

H = Hji + Hj2 + Heco + Heci + Hec2- (3) 

We recast the different parts of the Hamiltonian in Quantum Phase model. 
Hji = - Ej X J2i cos((f) i+1 - (pi), H J2 = - E J2 Ei cos(4> i+2 - <t>i). 

Hamiltonians Hji and Hj 2 are Josephson energy Hamiltonians respectively for nearest 
neighbor (NN) and next-nearest-neighbor (NNN) Josephson tunneling between the SQD. 

Heco = ^ Ei (-^^- - y) 2 , H E ci = E Z iY,i n i n-i+i, H E C2 = E Z2 J2ini n i+2 . 
H E co, H E ci, and H E C2 are respectively the Hamiltonians for on-site, NN and NNN charg- 
ing energies of the SQD. In the phase representation, (— igf:) is the operator representing 
the number of Cooper pairs at the ith dot, and thus it takes only the integer values (nj). 
Hamiltonian Heco accounts for the influence of gate voltage (eN ~ V g ). eN is the average 
dot charge induced by the gate voltage. When the ratio — * 0, the SQD array is in the 
insulating state having a gap of the width ~ Eco, since it costs an energy ~ Eco to change 
the number of pairs at any dot. The exception are the discrete points at iV = 2n + 1, where 
a dot with charge 2ne and 2(n + l)e has the same energy because gate charge compensates 
the charges of extra Cooper pair in the dot. On this degeneracy point, a small amount of 
Josephson coupling leads the system to the superconducting state. We are interested in 
analyzing the phases explicitly near this degeneracy point. 
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Now we want to recast our basic Hamiltonians in the spin language. During this process we 



follow Ref. 



H 



ji 



2E J1 J2 i (SjS i+1 -+h.c),H J2 = -2E J2 Ei(S^S i+2 - + h.c) 
Hecx — ^E Z \ ^2 S z S i+ i Z , 

i 

Hec2 = kEz2^TjSj S i+2 Z ■ 

i 

Here h = ( jV ~^ n ~ 1 ) allows tuning of the system to a degeneracy point by means of gate 
voltage. The phase diagram is periodic in iV with period 2. Here we only consider a 
single slab of the phase diagram, < N < 2. One can express spin chain systems to 
spinless fermions systems through the application of Jordan- Wigner transformation. In 
Jordan- Wigner transformation the relations between the spin and the electron creation and 

i-1 



annihilation operators are 

St = - 1/2 , S7 = (-l)V, exp[m J2 %] . $t = (S, 



t 

i J i 



J=-oo 



where rij = ipjipj is the fermion number at site j. Therefore, 



(4) 



i 

H J2 = -2E J2 £ (^ +2 ^ + h.c.)(4 +1 ^ +1 -l/2), 



(5) 
(6) 



H 



ECO 



-2hE C0 £ (^Vi-1/2), 

i 

Hecx = 4E Z1 - l/2)(^ +1 ^ +1 - 1/2), 

i 

H EC 2 = ±E Z2 (^-1/2)^2^2 -1/2). 



(7) 

(8) 
(9) 



In order to study the continuum field theory of these Hamiltonians, we recast the spinless 
fermions operators in terms of field operators by this relation 



ikpx 



(10) 
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where i/)r(x) and i^l{x) describe the second-quantized fields of right- and the left-moving 
fermions respectively, k F is Fermi wave vector. It reveals from Eq. 7 that the applied 
external gate voltage on the dot systems appears as a magnetic field in spin chain. So 
the different values of applied gate voltage promote the system in the different state of 
magnetization. In our system k F will depend on the applied gate voltage. The Fermi wave 
vector and the magnetization state (to) are related by the relation hp = |(1 — 2m). We 
would like to study the effect of gate voltage in our study so we keep k F as arbitrary. In our 
study we will mainly focus on the magnetization state with to = 0, 1/4, 1/2. We would like 
to express the fermionic fields in terms of bosonic field by the relation 

= JL= e"* " (11) 
V27ra 

r is denoting the chirality of the fermionic fields, right (1) or left movers (-1). The operators 
U r commutes with the bosonic field. U r of different species commute and U r of the same 
species anti-commute. field corresponds to the quantum fluctuations (bosonic) of spin and 
9 is the dual field of 0. They are related by the relations 4>r = 9 — <fi and <p L = 9 + 0. 
The model Hamiltonian after continuum field theoretical studies for arbitrary values of kp 
is 

H = H + AEz \ f : cos{A^fK(t){x) - (G - Ak F )x - 2k F a) : dx 
(2na) J 

+ AEz \ [ : cos(AVK(f)(x) +(G- Ak F )x - Ak F a) : dx (12) 
(2-ko) J 

2 



Ho = (f + ^-4£ J2 -^M [ dx [: (djf : + : (d^f : 

+ (16£ co - SE J2 - 4^-) / dx : (d x - d x cf>)(d x 9 + d x <f>) : (13) 
&co J 

H is the non-interacting part of the model Hamiltonian, v is the velocity of low energy 
excitations liquid parameter and the other is K (= /[ — 2 — - — J1 ^ 2 ]). 

Y E.n + %(4E Z1 - Tg^) 

2.1 CALCULATIONS AND RESULTS FOR COOPER PAIR STAIR CASE FOR 
SINGLE PAIR IN ALTERNATE SITES ( m = MAGNETIZATION PLATEAU) 



In this sub-section we discuss the Cooper pair stair case for a single pair in alternate sites. 
We have already mentioned that the SQD array system is nothing but the spin- 1/2 chain and 
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stair case of the SQD is nothing but the magnetization plateau. We will see that at m = 
(kp = it/2) there is an evidence of magnetization plateau. We would like to mention basic 
criteria for the appearance of magnetization plateau before we start our full swing discussion 
on stair case physics. The basic criteria for the appearance of plateau is the following: It is 
well known to us that the elementary excitation of an one dimensional spin system is gapless 
for half integer spin chain and gapped for integer spin chain. In the presence of magnetic field, 
it is possible for an integer spin chain to be gapless with partial magnetization and a half-odd- 
integer spin chain to show a gap above the ground state for appropriate values of the magnetic 



field. It has been shown by different groups [id, llll, ll2l, ll_3|, LUj that the magnetization 
of some systems can exhibit plateaus at certain nonzero values for some finite ranges of 
the magnetic field. The basic criteria for the appearance of magnetization plateau can be 
understood from the extension of Lieb-Schultz-Mattis theorem under a magnetic field. This 
implies that translationally invariant spin chains in an applied field can be gapful without 
breaking translation symmetry under the condition S — m = integer, where S and m are the 
spin and the magnetization state of the chain. In this gapful phase magnetization plateau 
occurs for quantized values of m. Fractional quantization can also occur, if accompanied by 
spontaneous breaking of translational symmetry. Fractional quantization can be understood 
from the S — m = non integer condition. In this situation system is either in the gapless 
low lying states or the degenerate ground state with spontaneous translational symmetry 
breaking in the thermodynamic limit. These conditions for the appearance of plateau are 
the necessary but not the sufficient condition. The nature and the occurrence of the plateau 

nn 

are also dependent on the nature of interaction present in the system [15L H6j. 

Our model Hamiltonian for m = magnetization plateau, in terms of bosonic fields, are 

H = H - A ^ Ezi ~ E 2 Z2) f cos{AVK(j){x) ) dx (14) 



(27ra)' 

In the spin chain, system is in the Neel-Ising state when XXZ anisotropy is greater than 
1 (4-Exi > 1). This state is doubly degenerate state and S — m = 1/2 , so this phase 
of the system satisfy the all criteria of fractionally quantized magnetization plateau. This 
is the relevant physics when the NNN exchange interaction is less than a critical value. 
When NNN exchange interaction exceed some critical value, ground state of the system is 



dimerized and doub 
of adjacent sites [17, 



y degenerate. The dimerized ground state is the product of spin singlet 



181 ] . In both ordering the lattice translational symmetry is breaking. 
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In our clean SQD system, we have predicted the CDW wave state with lattice translational 
symmetry breaking is alike to Neel-Ising state of one dimensional spin systems. This CDW 
state is the one-dimensional counter part of the two-dimensional checker board phase. One 
can find the dimer order density wave for the lower values of Cooper pair density, when one 
consider further long range interaction between the SQD. To get the dimer-order density 
in the presence of NNN interactions between dots, we shall have to consider the charge 
degeneracy point at higher values of half-integer Cooper pair density like n = 3/2. This 
Cooper pair stair case is originating due to the intersite Coulomb repulsion, so it is related 
with the Coulomb blocked phenomena of the system. 

2.2 CALCULATIONS AND RESULTS FOR COOPER PAIR STAIR CASE FOR 
SINGLE PAIR IN EVERY TWO SITES (m = 1/4 MAGNETIZATION PLATEAU) 

In this sub-section we discuss the occurrence of finite magnetization plateau (m = 1/4). 
This plateau state correspondence to the Cooper pair stair case for single Cooper pair in 
every two sites. For weak dimerization, i.e., when there is no reduction of unit cell size, 
reciprocal lattice vector is 2ir then the Hamiltonian of the system reduced to Eq. 15. None 
of the sine-Gordon coupling terms will provide the significant contributions for this plateau 
state, due to the oscillatory nature of the integrand. Effective Hamiltonian of the system 
for this state is the following: 



For strong dimerization state, reciprocal lattice vector G reduced to tx due to the reduction 
of the unit cell. The model Hamiltonian of the system reduced to 



where E Z12 = \J E Z \ + E Z2 2 and f3 = tan" 1 ^^. To study the renormalization group 
flow diagram, we shall have to construct the renormalization group (RG) equations. The 
RG equations of H 2 is 




(15) 




(16) 



dKjl) 
dl 



E zl2 2 K{i) 2 
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= (2-AK(l))E zl2 (17) 

In Fig. 1, we present the renormalization group flow diagram. RG flow lines, flowing off to 
E Z \2 — > oo are an Ising-Neel state, an anti-ferromagnetic ordering is setup in the z direction, 
i.e., the system is in the CDW wave state with only one boson in every two sites. This is 
the, Bl phase in the RG flow diagram. When the RG flow lines, flowing off to Ezu — oo 
are the dimer order density waves. In this state two Cooper pair of adjacent sites are bound 
with each other. This is the insulating state of the system. In spin chain, dimerized ground 
state is the analogous state of this phase. This phase is depicted as a B2 phase in the RG 
flow diagram. Phase D is the repulsive Luttinger liquid phase, i.e., the gap- less phase due to 
the irrelevancy of sine-Gordon coupling terms. Phase E is the superconducting phase of the 
system. This Cooper pair stair case is originating due to the intersite Coulomb repulsion, 
so it is related with the Coulomb blocked phenomena of the system. 



2.3 CALCULATIONS AND RESULTS FOR COOPER PAIR STAIR CASE AT 
EMPTY BAND LIMIT (m = \ MAGNETIZATION PLATEAU): 

Now we discuss the saturation plateau at m — \ (/cp = 0). K F = implies that 
the band is empty and the dispersion is not linear, so the validity of the continuum field 
theory is questionable. Values of the two Luttinger liquid parameters, v and K, are and 
1 respectively. It also implies that none of the sine-Gordon coupling terms become relevant 
in this parameter space. Saturation plateaus are only appearing due to very high values 
of magnetic field. In this plateau, system is in ferromagnetic ground state and restore the 
lattice translational symmetry. We think, this is the classical phase of the system. 

Absence of other fractionally quantized magnetization plateaus: Here we present the 
possible explanation for the absence of other fractionally quantized magnetization plateaus 
(like |, | etc): A careful examination of Eq. 12 reveals that to get a non oscillatory con- 
tribution from Hamiltonian one has to be satisfied = G condition but this condition is 
not fulfilled for these plateaus. The integrand of this sine-Gordon coupling terms contain 
an oscillatory factor that leads to a vanishing contribution. The other criteria is that non 
vanishing sine-Gordon coupling term should be relevant. 
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3. MODEL HAMILTON! AN AND PHYSICAL ANALYSIS OF COOPER PAIR 
STAIR CASE OF SITE DEPENDENT COUPLINGS, SUPERCONDUCTING 
QUANTUM DOTS 



At first we write down the model Hamiltonian of SQD system with site dependent Joseph- 
son couplings, on-site charging energies and intersite interactions in presence of gate voltage. 

H = Hji + H J2 + H EC o + H EC \ + H EC 2- (18) 

We recast the different parts of the Hamiltonian in Quantum Phase model. 

Hn = -Ej 1 '£ i (l-(-l) i 6 1 )cos((i>i +1 -<j)i),Hj2= -E J2 EiCos( y ^i + 2-(pi). 

Here NN Josephson couplings are different on alternate sites. Eji{\ + Si) and Eji(1 — Si) 
are the Josephson coupling strength for odd and even site respectively. 



Heco = ^ a E,(l-(-l) i 5 2 )(-^; 



d JV\2 

2 I ' 



Eco(l+S 2 ) and E C o(l— S 2 ) are the on site charging energies for odd and even site respectively. 

H E ci = EziYK 1 ~ (-l)%)ni n i+ i, 

i 

Ezi(l + S 3 ) and E Z i(l — S 3 ) are the inter-site charging energies for odd and even site re- 
spectively. All <f s are the deviational parameter from the clean superconducting quantum 
dots. 

H F C2 = E Z 2 ^ Uj n i+ 2, 

i 

Heco, Heci, and H E C2 are respectively the Hamiltonians for on-site, NN and NNN charging 
energies. Following the prescription of previous sections, one can write the continuum model 
Hamiltonian. 

H = H - 2E n J dx : sin(2^K (f){x) - (2k F - n)x) : 

+ ^ / : cos{AVK(/){x) - (G - Ak F )x - 2k F a) : dx 



(2na) 
| 4:E Z iS 3 

(2vra) 2 

| 2hE C0 

na 
AE 



J (-If : cos(A^fK(j)(x) - (G - Ak F )x - 2k F a) : dx 
J (-l) x : sin(2VK(j>(x) + 2k F x) : dx 



+ ^ / : cos{AVK<f){x) + (G - Ak F )x - 4k F a) : dx (19) 

(2ira) J 



3.1 CALCULATION AND RESULTS FOR COOPER PAIR STAIR CASE OF A 
SITE DEPENDENT SQD ARRAY FOR SINGLE COOPER PAIR IN ALTERNATE 
SITES (m = MAGNETIZATION PLATEAU) 

In this sub section, we find the evidence of Cooper pair stair case for single pair in 
alternate site for site dependent SQD system. The effective Hamiltonian reduce to 



At around the charge degeneracy point, the system is in the mixed state (M), i.e., the 
simultaneous presence of dimer density wave and staggered phase. Our model Hamiltonian 
consists of three sine-Gordon couplings. The first one arises due to site dependent variations 
of NN Josephson couplings, it yields the dimerized phase of the system. The anamolous 
scaling dimension of this term is 2K. This phase is spontaneous, i.e., infinitesimal variation 
of NN Josephson coupling is sufficient to produce this state. We have predicted in the 
previous sections of this work that the stair case is appearing due to the Coulomb blocked 
phenomena. Here we observe that the stair case may also occur for the Josephson tunneling. 
This prediction was absent in the previous studies. This observation of dimerized phase is in 
contrast with the clean superconducting quantum dots array, where it appears when NNN 
Josephson coupling exceed some critical value. The second sine-Gordon coupling term arises 
due to the site dependent applied gate voltage. The anamolous scaling dimension of this 
term is 2K. It yields the staggered phase of the system. The system is in the mixed phase 
(M) when both of the couplings are in equal magnitude otherwise the system is in any one of 
the state of this mixed phase, depending on the strength of couplings. We have ignored the 
3rd sine-Gordon coupling term of the Hamiltonian, because the anamolous scaling dimension 
is much larger than the other two. 




^4- [ dx : cos{A\f~K(b{x 

(2vra) 2 J 1 



(20) 
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3.2 CALCULATION AND RESULTS FOR COOPER PAIR STAIR CASE OF A SITE 
DEPENDENT SQD ARRAY FOR SINGLE COOPER PAIR IN EVERY TWO SITES 
(m = 1/4 MAGNETIZATION PLATEAU) AND EMPTY BAND LIMIT (m = 1/2 
MAGNETIZATION PLATEAU) 

In this sub section, we find the evidence of Cooper pair stair case for single pair in every 
two sites and also for empty band limit of site dependent SQD system. 

m — 1/4: in the weak dimerization limit, Hamiltonian reduced to 

H = H - 4Ezi t / dx : sin{WK<f>{x)) : (21) 

The sine-Gordon coupling term became relevant when K < 1/2. If we analyze the expression 
for Luttinger liquid parameter then we get the following quadratic equation of the system 
parameters to achieve this stair case phase for weak dimerization limit. 

512 

E n 2 - E n {A - —E zl ) - AE Z1 < 

This phase is in contrast with the stair case phase of the homogeneous SQD. In weak 
dimerization, there is no evidence of stair case phase for homogeneous SQD system. This 
phase of the system satisfies the all criteria of fractional quantization. Now we discuss for 
stronger strength of dimerization, i.e., when the dimerization strength exceed some critical 
value, at this point reciprocal lattice vector G reduced from 2n to tt. The Hamiltonian of 
the system reduced to 

H = H - AEzi2 2 f dx : cos(Ay/K(j>(x) + (3) : (22) 
{2-kol) J 

Expression for E Z n and (3 has given in section 2.2. This equation is identical to Eq. 
16. So the physics of the system is the same as Fig. 1. Hence the physical behavior of 
the system in strong dimerization limit is the same for homogeneous and site dependent SQD. 

m — 1/2 and other fractionally quantized plateaus: In these phases non of the sine- 
Gordon coupling terms are relevant. So there is no plateau phases for this states of the 
system. Hence the physical behavior of the system in strong dimerization limit is the same 
for homogeneous and site dependent SQD. 
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FIG. 1: Renormalization group flow diagram of the SQD array. We have depicted the different 
phases of the model Hamiltonians by: Bl. Charge-density wave (CDW), B2. Dimer-order density 
wave; D. Second kind of Repulsive Luttinger liquid; E. Superconducting phase. 

4. CONCLUSIONS 

We have predicted the Cooper pair stair case phenomena for homogeneous and also for 
the site dependent SQD system. We have predicted interesting physics for different stair 
case. We also conclude that the appearance of Cooper pair stair case is not only related with 
the Coulomb blocked phenomena but also related with the Josephson junction tunneling. 

The author would like to acknowledge the Center for Condensed Matter Theory of IISc 
for providing the working space. 
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